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Abstract 

We report the first measurement of the forward-backward asymmetry and the 
ratios of Wilson coefficients Ag/Aj and A\^jA-j in B — > K*£ + £~, obtained 
using 386M BB pairs that were collected at the T(4S) resonance with the Belle 
detector at the KEKB asymmetric-energy e + e~ collider. We also summarise 
the results obtained by Belle in searches for purely leptonic B decays, with 
emphasis on their impact on models of physics beyond the Standard Model. 



1 Introduction 



The Belle detector -0, operating at the KEKB e + e~ collider ^) at a centre-of- 
mass energy corresponding to the mass of the T(4S*) resonance, has accumu- 
lated as of May 2006 a data sample corresponding to more than 500M BB pairs. 
Such a large sample of B mesons allows the study of rare B decays, which are 
among the cleanest probes of the flavour sector of the Standard Model (SM) 
available to present experiments. In particular, these decays have the potential 
of revealing the existence of new particles and couplings not present in the SM, 
but predicted by several models of physics Beyond the SM (BSM). 

In the last few years Belle has reported several measurements that con- 
strain the parameter space of BSM theories. In particular, this review will 
concentrate on the measurement of forward-backward asymmetry and of ratios 
of Wilson coefficients in B — > K*£ + £~ described in Section 2, and on the 
searches for purely leptonic B decays of the type B — * I v, detailed in Section 3, 
and B° -> l+£- (Section 4). 



2 Measurement of forward- backward asymmetry and ratios of Wil- 
son coefficients in B — ► K*£ + £~ 

Flavor-changing neutral current b — > s processes are forbidden at tree level 
and can only proceed via loop diagrams in the SM. Loops are sensitive to 
new physics effects via insertion of heavy particles in the internal lines; if new 
heavy particles can contribute to the decays, their amplitudes will interfere with 
the SM amplitudes and thereby modify the decay rate as well as differential 
distributions. 

Such contributions would therefore change the Wilson coefficients ^) that 
parametrise the strength of the short distance interactions. The b — ► s£ + £~ 
amplitude is described by the effective Wilson coefficients C* s , Cg S and Cf^ , 
whose terms have been calculated up to next-to-next-to-leading order (NNLO) 

in the SM ^) . A measurement of the forward-backward asymmetry and of the 
differential decay rate g(q 2 , 9) = d 2 T/dq 2 d cos 9 as functions of q 2 and 9 for B — > 

K*£ + £~ constrains the relative signs and magnitudes of these coefficients ^' 
Here q 2 is the squared invariant mass of the dilepton system, and 9 is the angle 
between the momenta of the negatively (positively) charged lepton and the B 
(B) meson in the dilepton rest frame. The forward-backward asymmetry is 
defined as 

, , 2^ j\sgn(cos9)g{q 2 ,9)dcos9 

Afbw ) = 1 • (1) 

j\g(q 2 ,9)dco S 9 



2.1 Selection of B -> if events 



The analysis is based on a 357 fb 1 data sample containing 386M BB pairs; 
the B + — ► K + £ + £~ mode is also studied, since it is expected to have a very 
small forward-backward asymmetry even in the presence of new physics ^) and 
can therefore be used as null test for the „4fb measurement. 

The following final states are used to reconstruct B candidates: K* a £ + £~, 
K*+ £+£-, and K+£+£~, with subdecays K*° -► K+w', K*+ -» if°7r+ and 
ivT+TT , ifg 7T+7T-, and 7T° -> 77. Hereafter, K*°£+£~ and K*+£+£- are 
combined and called K*£ + £~. Two variables defined in the centre-of-mass 
(CM) frame are used to select B candidates: the beam-energy constrained mass 
M bc = \J ^bcam ~~ Pb and tne energy difference AE = E B - E hcam , where p B 
and i^s are the measured CM momentum and energy of the B candidate, and 
-Steam is the CM beam energy. 

The dominant background consists of BB events where both B mesons 
decay semileptonically. This background is suppressed using missing energy 
and cos 6* B , where 6* B is the angle between the flight direction of the B me- 
son and the beam axis in the CM frame. The resonant contributions due to 
charmonia (J/ip and i[)(2S)) are rejected using the dilcpton invariant mass. 
The signal box is defined as \M^ C — ms| < 8 MeV/c 2 for both lepton modes 
and -55 (-35) Me V < AE < 35McV for the electron (muon) mode. An un- 
binned maximum-likelihood fit to the M^ c distribution is performed in order 
to determine the signal yield. The fit includes several event categories, i.e. 
signal, three types of cross-feeds and four background sources. In the fit, all 
background fractions except the dilepton background are fixed while the signal 
fraction is allowed to float. 

The fit yields 113.6 ± 13.0 and 96.0 ± 12.0 signal events for the K*£+£- 
and K + £ + fr modes, respectively. Data and fit results are shown in Figure 1. 

2.2 Extraction of „4fb and Wilson coefficients 

The B — > K*£ + £~ candidates in the signal box are used to measure the nor- 
malised double differential decay width. For the evaluation of the Wilson co- 
efficients, the NNLO Wilson coefficients (7, of Ref. ^) are used. Since the full 
NNLO calculation only exists for q 2 /ml < 0.25, we adopt the so-called partial 

NNLO calculation ^) for q 2 /ml > 0.25. The higher order terms in the Cj are 
fixed to the SM values while the leading terms Ai, with the exception of Af, 
are allowed to float. Since the branching fraction measurement of B — > X s j 
is consistent with the prediction for \Aj\ within the SM, A7 is fixed at the 
SM value, —0.330, or the sign-flipped value, +0.330. The fit parameters are 
therefore A9/A7 and A w /At : the SM predictions for A 9 and Aiq are 4.069 and 

-4.213, respectively ^). 




Figure 1: M bc distributions for the (a) B -> K*l+(r and (b) B — ► K+£+£~ 
samples. The solid and dashed curves are the fit results for the total and back- 
ground contributions, respectively. 

To extract these ratios, an unbinned maximum likelihood fit is performed 
to the events in the signal box with a probability density function that includes 
the normalised double differential decay width. 

The g 2 -integrated asymmetry *4fb is measured by determining the yield 
in forward and backward regions from a fit to the M bc distribution. After 
correcting for efficiency, the result is: 

A F b{B K*£+r) = 0.50 ±0.15 ±0.02, 
Afb(B+ -> K+e+r) = 0.10 ±0.14 ±0.01, (2) 

where the first uncertainty is statistical and the second is systematic. A large 
integrated asymmetry is observed for K*£ + £~ with a significance of 3.4o\ The 
result for K + £ + £~ is consistent with zero as expected. 

The fit results of ratios of Wilson coefficients are summarised in Ta- 
ble 1. Figure 2 shows the fit results projected onto the background-subtracted 
forward-backward asymmetry distribution in bins of q 2 . 

The fit results are consistent with the SM values Ag/A? = —12.3 and 
Aio/Aj = 12.8. In Figure 3, we show confidence level (CL) contours in the 
(A9/A7, A10/A7) plane based on the fit likelihood smeared by the systematic 
uncertainty, which is assumed to have a Gaussian distribution. We also calcu- 
late an interval in AqAiq/A 2 - at the 95% CL for the allowed A? region, 

-14.0 x 10 2 < A 9 A W /A 2 7 < -26.4. (3) 



This result implies that the sign of A 9 A W must be negative, and the solutions 



Table 1: Ag/A 7 and A w /A 7 fit results for negative and positive A 7 values. The 
first uncertainty is statistical and the second is systematic. 





Negative A 7 


Positive A 7 


Ag/A 7 

A w /A 7 


-15.3±1;| ± 1.1 
10.3 ±^5 ± 1.8 


-16.3 +l\l± 1.4 
11.1 ±|;g ± 2.4 



in quadrants I and III of Figure 3 are excluded at 98.2% CL. Since solutions in 
both quadrants II and IV are allowed, we cannot determine the sign of A 7 A W . 

3 Search for B —> £ u 

Leptonic decays of charged B mesons proceed in the SM via the VF-mediated 
annihilation tree diagram, with a branching fraction given by: 

B(B+ - t*v e ) = ^P^mj (l ^f) f 2 B \V ub \*T B , (4) 
8tt V m B J 

where tb is the B meson lifetime, fs is the B decay constant and V u b an 
element of the CKM-matrix. These modes are thus very interesting because 
they give direct access to the product fs x V u b, from which one can extract a 
measurement of Jb- The SM expectation for B + — ► t + v t is around B(B + — > 
t + v t ) = 1.0 x 10~ ; decays to lighter leptons are helicity suppressed, by factors 
of 223 for muons and 10 7 for electrons. Allowing for decay amplitudes BSM, 
measurements of these processes give stringent limits on important parameters 
of such SM extensions, e.g. the mass of the charged Higgs boson and tan/3 (the 
ratio of vacuum expectation values of the two Higgs doublets) in the minimal 
supersymmetric SM (MSSM), or leptoquark masses in Pati-Salam models. 

3.1 B+ -> t+v t 

The search for B + — > t + v t y decays is based on the full reconstruction of one of 
the B mesons in the event, in the modes B + — > D^°h + and L)(*) L)i* )+ (h = 

n,K), selected in a sample of 275M BB events Decays of the types r — ► 
/z(e)i/P, ttv, 7T7r ^, and tttttti' are then searched for in the event remainders. The 

X A few weeks after this conference, the Belle collaboration has announced 
evidence for B + — > t + v t 

10) _ 

thus providing the first evidence of a purely 
leptonic B decay and the first direct determination of the decay constant 




Figure 2: Fit result for the negative A 7 solution (solid) projected onto the back- 
ground subtracted forward-backward asymmetry, and forward-backward asym- 
metry curves for several input parameters, including the effects of efficiency: 
A-j positive case (A7 = 0.330, Ag = 4.069, A10 = —4.213) (dashed), A±o pos- 
itive case (A 7 = -0.280, A 9 = 2.419, A w = 1.317) (dot-dashed) and both A 7 
and A w positive case (A 7 = 0.280, A 9 = 2.219, A w = 3.817J (dotted). The 
BSM scenarios shown by the dot-dashed and dotted curves are excluded by this 
measurement. 



final event selection is based on the extra energy present in the electromagnetic 
calorimeter. So-called double-tag events (i.e. fully reconstructed T(4S) decays) 
are used to validate the simulation of this quantity. The resulting upper limit 
on the branching fraction is B(B+ -> t+u t ) < 1.8 x 10~ 4 at 90% CL. In the 

two-Higgs doublet model, this limit can be translated 

12) 

in the exclusion of a 

region in the Higgs mass (m H + ) versus tan (5 plane. The Belle result is shown in 
Figure 4 together with constraints placed by other measurements ' 

14, 15). 



3.2 B+ -> n+Vp and B+ e+v e 

The electron and muon channels are studied in Belle by requiring one highly 
energetic lepton in the event and large missing energy and missing momentum 
corresponding to the undetected neutrino. All other particles in the event are 
used to reconstruct the companion B, which has to satisfy requirements on My, c 
and AE. The final selection variable used to isolate the signal is the lepton 
momentum in the B rest frame, defined by the momentum of the recoiling 
companion B. The latest Belle results are B(B+ -> ijl+v^) < 2.0 x 10~ 6 at 90% 

CL, based on a data sample of 140 ftT 1 16 ) and B(B+ e+v e ) < 5.4 x 10~ 6 

at 90% CL, based on 60 ftT 1 17 ). 
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Figure 3: Confidence level contours for negative A-j. Curves show la to 5a 
contours. The symbols show the fit (circle), SM (triangle), and A w -positive 
(star) cases. 

4 Search for B° —> £+ £~ 

The decay of B° to pairs of charged leptons can proceed in the SM via box 
diagrams or via annihilation into Z boson or photon. Similarly to their B^ 
counterparts, the predicted branching fractions are helicity suppressed for light 
leptons, down to about 1CT 10 for muons and 10~ 15 for electrons. Enhancements 
of these branching fractions are predicted by several BSM theories, such as high- 
tan (3 MSSM and some SUSY scenarios. The lepton flavour violating mode 
B° — > e ± [i T is forbidden in the SM and so any signal in this mode would be 
an undeniable sign of new physics. Belle has set the following 90% CL upper 
limits based on a data sample of 78 fb _1 18 ): B(B° -> ^ + /-0 < 1-6 x 1(T 7 , 
B{B° -> e+e") < 1.9xl(T 7 and B{B° -> e ± ^ T ) < 1.7xl(T 7 . The latter result 
can be interpreted as a lower bound on the Pati-Salam leptoquark mass 

19) : 

m LQ > 46 TcV/c 2 at 90% CL. 




Figure 4: The 90% C.L. exclusion boundaries in the [m#+ , tan /3] plane ob- 
tained from the observed upper limit on B(B + —>■ t + v), compared with other 
experimental searches. 
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